When scalar field is kinetically coupled to the Einstein tensor 
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We explore the cosmic evolution ol a scalar field with the kinetic term coupled to the Einstein 
tensor. We find that, in the absence of other matter sources or in the presence of only pressureless 
matter, the scalar behaves as pressureless matter and the sound speed of the scalar is vanishing. 
These properties enable the scalar field to be a candidate of cold dark matter. By also considering 
the scalar potential, we find the scalar field may play the role of both dark matter and dark energy. 
In this case, the equation of state of the scalar can cross the phantom divide, but this can lead 
to the sound speed becoming superluminal as it crosses the divide, and so is physically forbidden. 
Finally, if the kinetic term is coupled to more than one Einstein tensor, we find the equation of state 
is always approximately equal to -1 whether the potential is flat or not, and so the scalar may also 
be a candidate for the inflaton. 

PACS numbers: 98.80.Cq, 98.65.Dx 



I. INTRODUCTION 



Although at present there are no known fundamental 
scalar particles in nature, scalar fields play an important 
role in physics and cosmology. In physics, scalar fields are 
present in Jordan-Brans-Dicke theory as Jordan-Brans- 
Dicke scalar [l[; in Kaluza-Klein compactification the- 
ory as the radion Q, in the Standard Model of particle 
physics as the Higgs boson @ , in the low-energy limit of 
the superstring theory as the dilaton Q or tachyon [j| 
and so on. In cosmology, scalar fields are present as the 
inflaton Q to drive the inflation of the early Universe 
and currently as the quintessence or phantom [Toj 
fields, driving the acceleration the Universe. 

In general, the action of scalar-tensor theories of grav- 
ity are given by 



S = J d 4 x^y-[f(<f>,R, R^RT, R^XjR^) 

+K (<f>, d^<j>, V 2 0) + V (0) + S m ] , (1) 



where <j> is the scalar field, V{<j>) the scalar potential and 
R the Ricci scalar. Here / and K are arbitrary func- 
tions of the corresponding variables. These theories cover 
the f(R) modified gravity [llj, the Gauss-Bonnet grav- 
ity [l2[, the quintessence scalar, the phantom scalar, the 
quintom fields 13], the K-essence scalar theory the 
tachyon, dilaton, and so on. We note that these theories 
have already been investigated extensively. 

On the contrary, to our knowledge, the even more gen- 



eral scalar tensor theories of gravity such as 

s = J d 4 xV=g[f(<t>,R, R^BT, R^R^) 

+K (0, d^<t>, V 2 0, R^d^dvfr tRTdndvh ■ • •) 



+V (</>) + s„ 



(2) 



are has not been so well investigated. The new coupling 
between the derivative of scalar and the spacetime cur- 
vature may appear in some Kaluza-Klein theories flEh 
In 1993, Amendola [l8| studied the scalar-tensor 
theory with the Lagrangian linear in the Ricci scalar R, 
quadratic in 0, and containing terms as follows: 



Rd^c/) , R^(\)d v (f) , Rc/)W 2 
R ilv (f>d^d v (f> , d^Rd^cf) , V 2 i?0 . 



(3) 
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Amendola [HI investigated a cosmological model with 
the only derivative coupling term R fll ,d ,J " (j)d u (f> and pre- 
sented some analytical inflationary solutions. A general 
model containing Rd^d^cf) and R^yd^ <\>d v cj> has been 
discussed by Capozziello et al., O They showed 

that the de Sitter spacetime is an attractor solution in 
the model. In 2007, Daniel and Caldwell [2l| studied a 
theory with the derivative coupling term of R^d^^d 1, '</> 
and found the constraints on the coupling parameter with 
Solar system tests. 

In general, these scalar-tensor theories give both the 
Einstein equations and the equation of motion for the 
scalar in the form of fourth-order differential equations. 
However, recently Sushkov [22| showed that in the case 
of the kinetic term only coupled to the Einstein tensor, 
the equation of motion for the scalar is reduced to second 
order. Thus, from the point of view of physics, this theory 
can be interpreted as a "good" theory. The reason for 
this is very simple. It is well understood that tensors 
constructed from the metric tensor and its derivatives 
(up to second order), only the metric tenor and the 
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Einstein tensor are divergence free: 

V" = 0, G^ = 0. (4) 
Therefore for the Lagrangian 

2 = (f <T + f G"") d^d u ct> + V (0) , (5) 

with e and a constants, the equation of motion takes the 
form of 

(Eir + oGHW-^' 3 ". (6) 

which is a second order differential equation. 

In this paper, we will investigate the cosmic evolu- 
tion of a scalar field with the kinetic term coupled to 
more than one Einstein tensor. We find this scalar 
field presents us with two very interesting characteristics. 
When the kinetic term is coupled to only one Einstein 
tensor, and in the absence of any other matter sources 
or in the presence of only pressureless matter, the scalar 
behaves exactly as the pressureless matter. Thus it could 
be a candidate for cold dark matter. On the other hand, 
when the kinetic term is coupled to more than one Ein- 
stein tensor, the scalar field can have the equation of 
state to ~ — 1 over the whole history of the Universe. 
Therefore, the scalar field can play the role of a dynamic 
cosmological constant. 

The paper is organized as follows. In section II, we 
shall investigate the cosmic evolution of the scalar with 
the kinetic term coupled to one Einstein tensor. We find 
it to be a possible candidate for either cold dark matter 
or both cold dark matter and dark energy. In section III, 
we investigate the cosmic evolution of the scalar with the 
kinetic term coupled to more than one Einstein tensor. 
Section IV gives the conclusions and discussion. We shall 
use the system of units with G = c = h = k= l and the 
metric signature (— , +, +, +) throughout the paper. 

II. COUPLED TO ONE EINSTEIN TENSOR 

A. Equations of Motion 

The Lagrangian density of quintessence-phantom (or 
quintom) is given by 

2<& = \<r%W v * + V($) , (7) 

where is the metric tensor. Here e = +1 corresponds 
to quintessence and e = — 1 to a phantom field. Let's 
investigate the scalar field with the lagrangian density as 
follows 

^ = |g^<R0 + U(0) , (8) 

where 

= BP" - \g^R , (9) 



is the Einstein tensor, a is assumed to be a positive 
constant. We note that differently from quintessence or 
phantom fields, the kinetic term of this scalar field is cou- 
pled not to the metric tensor, but to the Einstein tensor. 
This modification leads to the scalar field behaving as 
the Einstein cosmological constant in Minkowski space- 
time due to the fact that G^ v = 0. Then from the action 




we obtain the Einstein equations 

G„ v = 87rT^ + 87rU 5 ^ + 87ra{-iV M 0V^i? 
+2V j8 0V ( „0i^ ) + V^V^i^ 
+V M V <T 0V„V <J 0- V A1 V I ^□0 
-±(V0) 2 G AlI ,+.g M^ ,[- \V a V^V a V p 4> 
+i(D0) 2 -V^V^^]} , (11) 

and the equation of motion for the scalar 

aG^dn&t - V = . (12) 

Here T^ v is the energy-momentum tensor for the matter 
fields derived from the action Sm. Prime denotes the 
derivative with respect to 4>. Due to the Bianchi identi- 
ties, 

G% = , (13) 

the equation of motion for the scalar field is equivalent 
to the conservation equation for other matter fields 

;i/ = 0. (14) 

We note that the equations of motion for a 
quintessence (or phantom) field is given by 

eg^d^cp - V = . (15) 

The gravitational field described the Solar system is the 
Schwarzschild solution which obeys 

GV = . (16) 

We conclude that as a test scalar, our Einstein tensor- 
coupled scalar field would behave as the Einstein cosmo- 
logical constant V = const in the Solar system, which is 
not the case for quintessence-phantom field. 

B. As Dark Matter 

In this subsection, we investigate whether the scalar 
can be a CDM candidate. Consider a spatially flat 
Fricdmann-Robertson- Walker (FRW) Universe 

ds 2 = -dt 2 + a 2 (dr 2 + r 2 dil 2 ) , (17) 
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where a is the scale factor. We model all other matter 
sources present in the Universe as perfect fluids. These 
matter sources can be baryonic matter, relativistic mat- 
ter and dark energy. We assume there is no interaction 
between the scalar field and other matter fields, other 
than by gravity. Then the Einstein equations are given 
by 



3H 2 = 



-- [2H + 3H' 



2H + 3H' = -8tt 



We define 



P 



p 



E 



Pi 



(18) 



(19) 



where pi, pi are the energy density and pressure of the 
i-th matter. 

The equation of motion for the scalar can be derived 
from above Einstein equations Eqs. (18) or from its equa- 
tion of motion Eq. (12) as follows 



Hcf> + (2B 



2H + 3H' 



= 



from which we obtain 



Co 



(20) 



(21) 



where cq is an integration constant. So the energy density 
of the scalar is 



9a 
P4> = ^T- 



r 2 



2 H 2 a e 



(22) 



The equation reveals that: in the stiff matter dominated 
Universe, namely H 2 oc 1/a 6 , the density of the scalar 
is nearly a constant; in the relativistic matter domi- 
nated Universe, namely H 2 oc 1/a 4 , the energy den- 
sity of the scalar scales as p$ oc 1/a 2 and the scalar 
has the equation of state —1/3 just like the curvature 
term; in the baryon matter dominated Universe, namely 
H 2 oc 1/a 3 , the scalar field behaves as pressureless mat- 
ter; in the cosmological constant dominated Universe, 
namely H 2 = const, the energy density of the scalar 
scales as p^ oc 1/a 6 and so has the equation of state +1. 
Now we have a Friedmann equation as follows 



3iT 



By setting 



V 2H 2 a e 



2 27ac 2 
^ 16tt 



P 



the Friedmann equation can be rewritten as 



3H Z = 8tt \p ■ 



4^2 



1 



(23) 



(24) 



(25) 



The energy density of the scalar turns out to be 




(26) 



We conclude that, in the absence of all other matter (p — 
0) or in the presence of baryonic matter (p oc 1/a 3 ), the 
scalar behaves exactly as pressureless dark matter with 
p<fr oc 1/a 3 . This enables us to propose this scalar field as 
a candidate for dark matter. Assuming that the energy 
density of the other matter fields is given by 



P 



PrO PbO 



A 



(27) 



where p r o, pbo are the present-day energy density of 
relativistic matter and baryon, respectively, and A is the 
Einstein cosmological constant. Then the pressure and 
the equation of state for the scalar take the form of 



P<p 



A 



PrO 
„4 



4£ 2 



(28) 



(29) 



We see again the scalar behaves exactly as pressureless 
matter in the presence of only baryon matter. Let pO 
denote the present-day energy density of the Universe. 
We define 



_ PrO_ 
PO 

Po ' 



n 



bO = 

pa 



PbO 

Po 



p_ 

po 



(30) 



Then the equation of state can be written in the dimen- 
sionless form 



1 



w 



Q 2 



(31) 



Observations constrain il^o = 0.04 , tt r o = 8.1 • 
10~ 5 , n x = 0.75 . We find k = 0.458 correspond- 
ing to the density fraction of scalar p<f>o/pO = 0-21 in the 
present-day Universe (23[. 

In Fig. [1] we plot the equation of state for the scalar. 
It shows that the scalar behaves as a curvature term 
(w = —1/3), pressureless matter (w = 0) and stiff mat- 
ter (w — 1) in the relativistic matter, pressureless mat- 
ter and cosmological constant dominated epoch, respec- 
tively. In Fig. [51 we plot the density fraction of the scalar 
during the evolution of the Universe. It is found that the 
scalar makes significant contribution only at the interval 
of In a = — 9 to In a = 0. This is not surprising since the 
scalar only tracks the pressureless matter. The circled 
line is for ACDM model. The difference between the two 
models is observationally allowed as shown in Fig. [3J 
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FIG. 1: The equation of state for the scalar (solid line) and 
the cold dark matter (circled line). The scalar behaves as a 
curvature term (w = —1/3), pressureless matter (w = 0) and 
stiff matter (to = 1) in radiation, scalar matter and cosmo- 
logical constant dominated epoch, respectively. 
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FIG. 2: The evolution of density fractions 



p<p+p 



of the scalar 



field (solid line) and pd/p of the cold dark matter (circled line) 
in the ACDM model. The scalar contributes significantly only 
at the interval of In a — —9 to In a = 0. 




FIG. 3: The evolution of energy density logp^ of the scalar 
field (solid line) and log pd of the cold dark matter (circled 
line) in the ACDM model. The difference is observationally 
allowed. 



C. As Dark Matter and Dark Energy 

In the previous subsection, we have shown the scalar 
can play the role of dark matter. In this subsection, 
we study whether the scalar can play the role of both 
dark matter and dark energy by taking account the scalar 
potential. Since a is assumed to be a positive constant, 
we may rescale <f> such that the Lagrangian density can 
be rewritten as 



2 = -G» v d ll <$>d v 4> + v(4>) 



(32) 



Then the Einstein equations and the equation of motion 
take the form of 



3H 2 = 8tt Qff 2 ^ 2 + V + Pb + p r ^j , 

l -(2H + 'iH 2 ) </> 2 + 



2H + 3H 2 



-V+ P r} , 



H 2 j> + H[2H + 3H 2 )6 + W 



-8?r 



2 -U 

3 H v 



(2H + 3H 2 ^j j> 







(33) 



where pj,, p r are the energy density of baryon matter 
and relativistic matter. p r is the pressure of relativistic 
matter. Now we investigate the evolution of the scalar 
field. We introduce the following dimensionless quantities 



/47T ; 

T 

V 



v = 



8irV 



8nVH 



3H 2 ' 
V"V 

= y>2 



_ l &nPb 

Z ~ V3iP ' 
, TV ee In a. (34) 



Here x 2 and y 2 represent the density parameters of the 
kinetic and potential terms respectively. Then the above 
equations can be written in the following autonomous 
form 



dx 
dN 

dy 
dN 

dz 
dN 



\f& o X [ r- 2 

~^-y u + 1 — I — 2 -2v 6xy u 
2 1 + x z L 

-3y 2 + (l-x 2 -y 2 -z 2 )] , 

V6 3 y 

-— X yu+-y+ WT - r) 

-3y 2 + {l-x 2 -y 2 -z 2 )] , 

-2V6xy 2 u - 3y 2 



-2V6xy 2 



2(l + x 2 ) 



+ (l-x 2 -y 2 -z 2 )] 



= — VETxu 2 + VExu 2 

dN 



2(l + x 2 ) 

-3y 2 + (l-x 2 -y 2 - 
together with a constraint equation 



-2V6xy 2 



\ -i 3u 



2 i 2 i 2 

x z +y z + z A 



8np r 
3H 2 



= 1. 



(35) 



(36) 
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FIG. 4: The phase plane for the exponential potential. The 
point (0, 0) corresponds to the radiation dominated epoch and 
the circled arc (x 2 + y 2 = 1) corresponds to the scalar domi- 
nated epoch. The point (0, 0) is unstable and the point (0 ,1) 
is stable and an attractor. These trajectories show that the 
Universe always evolves from the radiation dominated epoch 
to the scalar field dominated epoch and ends at the scalar 
potential dominated epoch. 



The equation of state w and the fraction of the energy 
density il(f) for the scalar field are 



fid, = 



V<t> 



{l + x 2 ){Zx 2 + Zy 2 ) 



2V6xy 2 u + 3y 2 



(l-x 2 -y 2 -z 2 )] - 



2y/Qxy 2 u + 3y 2 
3a; 2 + 3y 2 



2 , 2 
Jh 2 ~ = X +V ■ 



(37) 



As an example, we have considered the exponential po- 
tential V = with £ a positive constant. Physi- 
cally, the scalar would roll down the potential. There- 
fore we have x > 0, y > 0, z > 0, u > 0. We put 
il M = 0.04 , fJ r0 = 8.1 • 10~ 5 , n x = 0.75 . 

In Table HI we present the properties of the critical 
points for the exponential potential. The point (a) cor- 
responds to the relativistic matter dominated epoch and 
the point is unstable. In this epoch, the scalar has the 
equation of state —1/3. The line (b) corresponds to the 
scalar (kinetic energy)-plus-baryon dominated epoch and 
it is a saddle line. In this epoch, the scalar field behaves 
as dust with equation of state w — 0. Point (c) corre- 
sponds to the potential energy dominated epoch. It is 
stable and thus an attractor. In the epoch of (c), the 
expansion of the universe accelerates. 

In Fig. HI we plot the phase plane for the scalar with 
a range of different initial conditions. The point (0, 0) 
corresponds to the radiation dominated epoch and the 



circled arc (x 



y 



1) corresponds to the scalar domi- 



nated epoch. The point (0, 0) is unstable and the point 
(0 ,1) is stable and an attractor. These trajectories show 
that the Universe always evolves from the radiation domi- 
nated epoch to the scalar field dominated epoch and ends 
at the scalar potential dominated epoch. 

In Fig. [SI we plot the density fraction of scalar field 



FIG. 5: The solid line represents the density fraction of the 
scalar. The circled line represents the density fraction for 
the mixture of cold dark matter and cosmological constant in 
ACDM model. 




FIG. 6: The solid line represents the evolution of the equation 
of state for the scalar. The circled line represents the equation 
of state for the mixture of cold dark matter and cosmological 
constant in ACDM model. 



and the density faction of dark matter and cosmological 
constant in ACDM model. The solid line represents the 
density fraction of the scalar. The circled line represents 
the density fraction for the mixture of cold dark matter 
and cosmological constant in ACDM model. In Fig. [51 
we plot the equation of state for the scalar and for the 
mixture of cold dark matter and cosmological constant in 
ACDM model. The solid line represents the evolution of 
the equation of state for the scalar. The circled line rep- 
resents the equation of state for the mixture of cold dark 
matter and cosmological constant in ACDM model. In 
both cases, we set fijo = 0-04 , O r o = 8.1 • 10~ 5 , VL\ — 
0.75 . 



D. Crossing the Phantom Divide 

In the absence of baryon and relativistic matter, we 
find the scalar can cross the phantom divide (w = — 1). 
To demonstrate this, we write down the autonomous sys- 
tem of equations 
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Name 


X 


V 


z 


M 


Existence 


Stability 




w 


(a) 














All C, 


Unstable node 





-1/3 


(b) 







z 





All C, 


Saddle line 







(c) 





1 








All C 


Stable node 


1 


-1 



TABLE I: The properties of the critical points for the scalar for the exponential potential. 



dx v6 / ,x x 

= 1-r «H ;r 

dN 2 V ' 1 + x 2 



0.25 i 



-3(1-^)] , 



aiv 2(l + aH) 



-2V6x (1 - 
[-3(1 



a; 2 ) u 



-2V6X (1 



x 2 ) « 



3// 

T 



together with the constraint equation 

x 2 + y 2 = 1. 
The equation of state w is given by 



w 



P<i> 



1 2V/6U2; 3 - 2yJ§xy - 3 + 3x 2 
3 1 + x 2 



(38) 
(39) 

(40) 



In Table [HI we present the properties of the critical 
points for the exponential potential. The points (a) cor- 
responds to the kinetic energy dominated epoch and this 
point is unstable. In this epoch, the scalar has the equa- 
tion of state 0. The points (b) corresponds to the epoch 
where both the kinetic energy and potential energy van- 
ishes, and this point is stable and an attractor. In this 
epoch, the scalar field behaves as the cosmological con- 
stant. 

In Fig. [7] and Fig. [FJ we plot the phase plane and the 
behavior of the equation of state during phantom cross- 
ing. Although it can cross the phantom divide, in the 
next subsection, we show that the corresponding sound 
speed of the scalar is either imaginary or superluminal. 
So this phantom crossing behavior may be physically for- 
bidden. 



E. Sound Speed of the Scalar 

The scalar field is constant in Minkowski spacetime due 
to = 0. So it is trivial to consider the perturbation 
theory of the field in Minkowski spacetime. Here we con- 
sider the field perturbation in the background of FRW 
Universe. In this paper, we have assumed the scalar is 
not coupled to the Standard Model particles but only to 
gravity. Therefore, in order to derive the equation of mo- 
tion for the scalar field perturbation, we start from the 
Lagrangian density: 



,<t> + v(<f>) 



(41) 




FIG. 7: The phase plane for the exponential potential. The 
point (1, 0) corresponds to the kinetic energy dominated 
epoch and it is unstable. The point (0 ,0) corresponds to the 
potential energy dominated epoch and it is stable. These tra- 
jectories show that the scalar always evolves from the kinetic 
energy dominated epoch to the potential energy dominated 
epoch. 
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FIG. 8: The equation of state for the scalar with two dif- 
ferent initial conditions a:(0) = 0.1, u(0) = 0.1 and x(0) = 
0.1, w(0) = 0.9. For the condition x(0) = 0.1, «(0) = 0.9, the 
scalar can across the phantom divide. 



It is convenient to work in Newtonian gauge. In the ab- 
sence of anisotropic stress, and for scalar perturbations 
only, the perturbed FRW metric can be written in the 
form 



ds 2 



(1 + 2$) dt 2 + a (t) 2 (1 - 2$) (dr 2 + r 2 dfl 2 ) (,42) 



where $ is the gauge invariant Newtonian potential. The 
potential characterizes the metric perturbations. 
We proceed to linear order by perturbing the field 



(43) 
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Name 


X 


it 


Existence 


Stability | f2^ | w 


(a) 


1 





All C 


Unstable node | Kinetic Energy Dominated | 


(b) 








All C 


Stable node | scalar vanishes |-1 



TABLE II: The properties of the critical points for the scalar for the exponential potential. 



We find the explicit form of the perturbation equation 
has been given in Ref. [2J|: 



1 



a QU -s — [7 = 0, 



(44) 



where A is the three dimensional Laplace operator, s is 
the sound speed squared and 



[/ = $ 



H 



Q 



2F+Q b 



Qe 



Qa 



Q d + ^i 



1 + 



2F+Q b 



( 2F+Q b ) Of 



2F+Q b 



F E 








Qa 


2 




= - V , 


Qb 


9^ ' 


Qc 


EE -H 2 ^ 2 




3 ^ 


Qd 


2 • - 9 

9 v ' 




4 . /.. 


Qe 




Qs 


- - V . 

9^ 



(45) 

Using the variables defined in the subsection II-D, we find 
the square of the sound speed can be written as 



1 r 

——^ — ~s , n — 9w — 30 
3 (x 2 + 1) {x 2 - 3) L 

+4x (x 2 - l) (x + 3V6u - V6 



wx~ + 9wx 

2^ 



ux 



(46) 



Then using the autonomous system of equations Eq. 
we can plot the evolution of the square of the sound 
speed. 

In Fig. IH1 we plot the behavior of the sound 
speed squared. It shows that, in the kinetic term 
(G fJ ,uV fl <t>V / <t>) dominated Universe, the sound speed of 
scalar always vanishes. Then the vanishing of both 
pressure and sound speed enable the scalar (without 
a scalar potential) to be a plausible candidate of cold 




-0.5- 



FIG. 9: Evolution of the sound speed squared s of the scalar 
with two different initial conditions x(0) = 0.1, u(0) = 0.1 and 
x(0) = 0.1, u(0) = 0.9. For the condition x(0) = 0.1, u(0) = 
0.9, the sound speed yfs of the scalar is either imaginary or 
superluminal in the vicinity of phantom divide crossing (see 
Fig. [8}. Therefore, the phantom divide crossing behavior may 
be physically forbidden. 



dark matter. For some initial conditions, for example, 
x(0) — 0.1, w(0) = 0.9, the sound speed would be ei- 
ther imaginary or superluminal in the vicinity of phan- 
tom divide crossing (see Fig. [5]). So this behavior may 
be physically forbidden. But for some initial conditions, 
for example, x(0) = 0.1, u(Q) = 0.1, the sound speed 
of the scalar is always smaller than the speed of light. 
Therefore, it is viable that the scalar field (with a scalar 
potential term) is a potential candidate for dark matter 
and dark energy. 



III. COUPLED TO n EINSTEIN TENSORS 

A. Without Potential 

In section II, we investigated the cosmic evolution of 
the scalar field coupled to a single Einstein tensor. In 
this section, we shall extend this coupling to n Einstein 
tensors. Consider the following Lagrangian density of the 
scalar field, which is coupled to gravitation 

& = ^ ai GZG% ■ ■ ■ GlizlG^d^d^ , (47) 
where 



ff" 1 = - -g^R 

2 y 



(48) 



is the Einstein tensor. A is assumed to be a positive con- 
stant. By considering variation of the Lagrangian density 
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with respect to the field <f>, we obtain the following equa- 
tion of motion 



V, 



gzi ■ ■ •G""- 2 G"- 1,y a> =0 



(49) 



In general, it is a third order differential equation. But in 
the background of Friedmann- Robertson- Walker (FRW) 
Universe, it reduces to a second order differential equa- 
tion. We also model all other matter sources present in 
the Universe as perfect fluids. In this case the Einstein 
equations are given by 



3H 2 
2H + 3H 2 



8tt (p^ + p) , 
-8?r (p + p) 



(50) 



where the energy density and pressure of the scalar field 
are given by 



= \(2n+l)(3H 2 ) n <j> 2 , 

= \(2n + l)(3H 2 ) n <j> 2 ( 1 + 



2nH 
3H 2 



(51) 



p, p are defined by Eq. (19). 

The equation of motion of the scalar field derived from 
Eq. ((49]) takes the form 



(52) 



H(j> + 2nH + 3H Z U = 



This equation of motion is equivalent to the energy con- 
servation equation of the scalar. From the equation of 
motion we obtain 



Co 



H 2r, 



(53) 



with cO an integration constant. This gives us the follow- 
ing Friedmann equation, 



3H 1 



8tt 



Ac§3" (2n + l) 



(54) 



So the energy density of the scalar field scales as 
H~ 2n a~ 6 . In the radiation dominated epoch, we have 
H 2 ~ a~ 4 , and so the energy density of scalar field would 
scale as a 4n ~ 6 . It behaves as a phantom field, growing 
in density with the expansion, when n > 2. In the mat- 
ter dominated epoch H 2 ~ a~ 3 , so the energy density of 
scalar field scales as a 3 ™~ 6 . It behaves as cosmological 
constant for n = 2 and as a phantom field for n > 3. In 
total, when n > 3 the scalar behaves as phantom field in 
over the whole history of the Universe. 

Taking into account pressureless matter (dark matter 
plus baryonic matter) and relativistic matter, we obtain 
the Friedmann equation in the following dimensionless 
form 

h 2 = Cl r0 a- 4 + n m0 a- 3 + n X oa~ 6 h- 2n , (55) 
In Fig. 1101 we plot the deceleration parameter q: 



(56) 





FIG. 10: The deceleration parameter for the scalar (solid line) 
model and the ACDM (circled line) model. For the scalar 
model, from left to right, we set n = 6, 9, 12, 15, 20, 30, 
respectively. 



-0.7 




0.4 z -0.2 



FIG. 11: The equation of the state of the scalar field 
as a function of redshift. From top down, we set n = 
6, 9, 12, 15, 20, 30, respectively. They smoothly cross the 
phantom divide. 



as a function redshift for different values of n. We find 
the scalar model can really enhance the acceleration of 
the Universe. 

In Fig. [11] and Fig. [12l we plot the equation of state 
for the scalar. In Fig. [T3l we plot the energy den- 
sity flxo a ~ 6 H~ 2n of the scalar and the energy density 



■.oa 



of matter, respectively. For different values of 



n, the energy density is always asymptotically approach- 
ing zero at higher redshifts. For n > 9, the density is 
vanishing at redshifts greater than one. Thus by simply 
choosing a relative large n we obtain a dark energy neg- 
ligible at redshifts greater than one without fine tuning 
of initial conditions. We know the energy density of the 
usual phantom field grows during most of the history of 
the Universe, so adjusting it to become important very 
recently requires even more fine tuning of the initial con- 
ditions. However, in above models, for a very large range 
values of n [n > 9), the scalar field becomes important 
very recently, as shown in Fig. 1121 
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z 1 



FIG. 12: The equation of the state of the scalar field 
as a function of redshift. From top down, we set n — 
6, 9, 12, 15, 20, 30, respectively. They asymptotically ap- 
proach the constant equation of state at higher redshifts. 




FIG. 13: The energy density Qxoa-~ 6 h~ 2n of the scalar field 
(solid line) and the energy density ilmOa -3 of matter (circled 
line). The crossed line is for the cosmological constant. From 
right to left, we set n = 6, 9, 12, 15, 20, 30, respectively. 
They all asymptotically approach zero at higher redshifts. For 
n > 9, the density of the scalar is vanishing at redshifts greater 
than unit one. 



B. With Potential 



In this section, we investigate the cosmic behavior of 
the scalar field in the presence of a scalar potential. The 
Lagrangian density is given by 



* ~ 2(2n+l)3« G " lG ^ ' ' ' Gan ~ lUd ^ 

+v{4>) , 



(57) 



where n a positive integer. Then the Einstein equations 
are given by 



2H + 3H 2 



8?r (p$ + p) , 
-8tt (p + p) 



(58) 



where the energy density and pressure of the scalar field 
are given by 



= lu 2n 6 2 + V 



P0 



-H 



2nl2 



2nil\ 
3H 2 J 



2n 



V'^H' 1 - V .(59) 



The equation of motion of the scalar field derived from 
Eq. (gSJ takes the form 



3H<j> 1 



Ml 



2nH 



2n + l 



V =0 



V 3H 2 J H 2 
We introduce the following dimensionless quantities 



-H r 



V 



inVH r 



y 

r : 



8ttV_ 
3H 2 

v"v 



y>2 



(60) 



3H 2 ' 
N = In a. (61) 



Then the above equations can be written in the following 
autonomous form 



dx 

In 



dy 
dN 



dz 

In 

du 
dN 



v6 , „ o (n + 1) x r o 

_(2n + l)A- 3a : + A_^_ [4 + 2 ^ 

— Ay 2 — z 2 — 2\f§nxy 2 u 

-—xyu + — — ^ [4 + 2x 2 

2 y 2(1 + nx 2 ) 1 



-Ay 2 - z 2 - 2V6 



nxy 2 u 



-z H ; — rr \A + 2a; 2 

2 2(1 +w 2 ) L 



-Ay 2 - z 2 - 2V6nxy 2 
— V6Txu 2 + Vfixu 2 1 
-Ay 2 - z 2 - 2V6nxy 2 



n U: 



2(1 + nx 2 



[4 + 2x 2 



together with a constraint equation 
2 , 2 , 2 , 87rp r 

* +y + z +^2- = 1 - 

The equation of state w and the fraction of the energy 
density £]L for the scalar field are 



(62) 



(63) 



a, 



P0 (1 + nx 2 ) (3x 2 + 3y 2 ) 
+nx 2 — An + Any 2 + nz 2 ] — 



2 i 2 



2\/6n 2 xy 2 u + 3 

2\/6nxy 2 -it + 3y 2 
3x 2 + 3y 2 



3H 2 



1 + 

A + 2x 2 - Ay 2 - z 2 - 2^/Qnuxy 2 
~ 1+ 2{l + nx 2 ) ( ' 
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As an example, we have considered the exponential po- 
tential V — with £ a positive constant. 

In Table IIII| we present the properties of the critical 
points for the exponential potential. The point (a) cor- 
responds to the relativistic matter dominated epoch and 
this point is unstable. In this epoch, the scalar has the 
equation of state w = 1 — 4p . The point (b) corresponds 
to the kinetic energy of the scalar dominated epoch and 
is a saddle point. In this epoch, the scalar field has the 
equation of state w = fe^. The point (c) corresponds to 
the potential energy dominated epoch. It is stable and 
an attractor. In this epoch, the scalar has the equation of 
state w = — 1. The point (d) corresponds to the matter 
dominated epoch. It is a saddle point. In the epoch of 
(d), the scalar has the equation of state 1 — n. 

We note that when n = 1, the scalar behaves as: a 
curvature term in radiation epoch; pressureless matter in 
both the kinetic energy dominated epoch and the matter 
dominated epoch; a cosmological constant in potential 
energy dominated epoch. When n = 2, the equation of 
state is smaller than —1/3. 

In Fig. [14l we plot the phase portraits for the scalar 
with vast initial conditions. The point (0, 0) corresponds 
to the radiation dominated epoch and the circled arc 
(x 2 +y 2 = 1) corresponds to the scalar dominated epoch. 
The point (0, 0) is unstable and the point (0 ,1) is stable 
and an attractor. These trajectories show that the Uni- 
verse always evolves from the radiation dominated epoch 
to the scalar field dominated epoch and ends at the scalar 
potential dominated epoch. 

In Fig. [15l we plot the equation of state for the scalar 
when n — 2. As shown in Table \UJ\ the equation of state 
is smaller than —1/3. Buy choosing some smaller (, we 
have the equation of state w ~ — 1 in the total history 
of the Universe. In Fig. [TH1 we plot the deceleration of 
the Universe for the scalar model and ACDM model. We 
find the two models predict nearly the same behavior of 
the Universe from deceleration to acceleration. This is 
because the equation of state for dark energy is w ~ — 1 
over the entire history of the Universe, and so the energy 
density of this dark energy is nearly a constant. 




1.2 



FIG. 14: The phase plane for the exponential potential with 
n = 2. The point (0, 0) corresponds to the radiation domi- 
nated epoch and the circled arc (x 2 + y 2 — 1) corresponds to 
the scalar dominated epoch. The point (0, 0) is unstable and 
the point (0 ,1) is stable and thus an attractor. These tra- 
jectories show that the Universe always evolves from the ra- 
diation dominated epoch to the scalar field dominated epoch 
and ends at the scalar potential dominated epoch. 



-0.9- 




-15 



-10 



FIG. 15: Evolution of the equation of state for the scalar. 
The scalar can cross the phantom divide. But the equation 
of state w — — 1 over most of the history of the Universe. 



C. New Infiaton field 

In this section, we investigate the cosmic behavior of 
the scalar field in the absence of other matter sources. 
In this case, we have the following autonomous system of 
equations 



dx 
dN 



du 
dN 



V6 



(2n- 



-4(1-; 
-Ail-: 



1) (l - x 2 ) it - 3 
• 2 ) - 2\f&nx (1 - 



x 2 ) u 



\/6: 



n u 



xu 



2(1 



(n + l)x 
2(l + nx 2 ) 

[4 + 2a; 2 



[4 + 2x 2 



0.2- 



o- 



q 

-0.2 H 



-0.4 



-0.6 




2 v / 6nx (l 



x 2 ) u 



(65) 



FIG. 16: The evolution of decelerating parameters for ACDM 
model (pointed line) and the scalar model (solid line). Both 
models predict the transition redshift of the Universe from 
deceleration to acceleration at zT ~ 0.8. 
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Name 


X 


V 


z 


u 


Existence 


Stability 


^4> 


w 


(a) 














All C, 


Unstable node 





—I L 


(b) 


1 











All C, 


Saddle line 


1 


1 — n 

l+n 


(c) 





1 








All C, 


Stable node 


1 


-1 


(d) 








1 





All C 


Saddle node 





1 - n 



TABLE III: The properties of the critical points for the scalar for the exponential potential. 



-0.98 




FIG. 17: The equation of state for n = 100. 



together with a constraint equation 
n 8ttV , 



3H 2 



(66) 



The equation of state w is 



3(1 + nx 2 ) 



2V&n 2 x(l- x 2 )u + 3 



-3na; 2 ] ^—nx (l — x 2 ) u — 1 , 



(67) 



As an example, we consider the exponential potential 
V = with ( a positive constant. 

In Table IIV1 we present the properties of the critical 
points for the exponential potential. The points (a) cor- 
responds to the kinetic energy dominated epoch and this 
point is unstable. In this epoch, the scalar has the equa- 
tion of state jt^- The points (b) corresponds to the 
scalar vanishes and this point is stable and an attractor. 
In this epoch, the scalar field also has the equation of 
state —1. So by choosing a large n we find the equation 
of state is always approximately equal to —1, without fine 
tuning the initial conditions. For the infiaton field, one 
generally assumes that the scalar potential is sufficiently 
flat to make the equation of state w ~ — 1. In this case, 
our scalar field can easily achieve this equation of state 
without the flatness condition. Therefore this scalar field 
could act as a new candidate for the infiaton. 

In Fig. [TT1 we plot the equation of state for n — 100. 



IV. CONCLUSION AND DISCUSSION 

In this paper, we have investigated the cosmic evolu- 
tion of the scalar field with the kinetic term coupled to 
n Einstein tensors. When n = 1, the equation of mo- 
tion is a second order differential equation, and so is a 
good theory. We find that the scalar has some interest- 
ing properties. Firstly, in the absence of other matter 
sources or in the presence of only pressureless matter, 
the scalar (without a scalar potential) behaves exactly as 
pressureless matter. Further more, the sound speed of 
the scalar perturbations is exactly zero in the structure 
formation era. These properties enable the scalar to be 
a potential candidate for cold dark matter. Secondly, in 
the presence of a scalar potential, the scalar field has the 
equation of state — 1 < w < 0. So the scalar may play 
the role of both cold dark matter and dark energy. It is 
found the sound speed is always smaller than the speed 
of light, and so it is physically viable. By investigating 
the dynamics of the scalar field in the absence of other 
matter sources, we find it can cross the phantom divide. 
But it may be physically forbidden due to the presence 
of a superluminal sound speed. Finally, if the kinetic 
term is coupled to more than one Einstein tensors, the 
equation of state is always approximately equals to — 1 
whether the potential is flat or not. Thus the scalar may 
also be a potential candidate for the infiaton field. 
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